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ABSTRACT 

We construct a numerical model of emission from minijets, localized flows driven by magnetic 
reconnection inside Poynting-fiux-dominated jets proposed to explain the ultrafast variability 
of blazars. The geometrical structure of the model consists of two wedge-like regions of rela- 
tivistically flowing gas, separated by a stationary shock. The dynamics is based on solutions 
of relativistic magnetic reconnection with a guide field from Lyubarsky (2005). Electron dis- 
tributions in each region are chosen to the match the pressure and density of the local plasma. 
Synchrotron emission from both regions is used to calculate Compton scattering, Compton 
drag and photon-photon opacity effects, with exact treatment of anisotropy and the Klein- 
Nishina regime. Radiative effects on plasma are taken into account, including the dependence 
of pressure on electron radiative losses and adiabatic heating of the flow decelerating under 
Compton drag. The results are applied to the July 2006 flare in the BL Lac object PKS 2155- 
304, with the aim of matching TeV flux measurements by H.E.S.S. with models that satisfy the 
variability constraints, while keeping X-ray emission below simultaneous Chandra observa- 
tions. We find that models of isolated minijets with a significant guide field overproduce X-ray 
emission, and that we must take into account the radiative interaction of oppositely-oriented 
minijets in order to achieve a high enough dominance by Comptonized TeV radiation. We 
argue that such interactions are likely to occur in a jet where there is substantial internal re- 
connection, producing a large number of misaligned minijets. Finally, we show that large jet 
magnetizations are indeed required to satisfy all observational constraints and that the effec- 
tive Lorentz factor of the minijet plasma has to be larger than 50, in agreement with earlier 
one-zone estimates. 

Key words: magnetic reconnection - radiation mechanisms: non-thermal - galaxies: active 
- galaxies: jets - BL Lacertae objects: individual: PKS 2155-304 



1 INTRODUCTION 

Relativistic jets produced in Active Galactic Nuclei (AGNs) are 
thought to be launched as cold and highly-magnetized outflows, 
the energetics of which is dominated by Poynting flux. Under such 
conditions shocks are very inefficient, however energy dissipation 
can be provided by relativistic magnetic reconnection. This idea 
led to the development of the minijets model, which has been pro- 
posed to explain extre mely fast variability o bserved in TeV blazar s 
JGiannios et alj2009T) and radio galaxy M87 JGiannios et al.l201fj) . 
Minijets are perpendicular relativistic flows within a relativistic jet 
and as such attain very high Lorentz factors V > 50 required 
to cir cumvent the gamma-ray opacity problem l IBegelman et alj 
120081) . A number of alternative solutions for this 'Lo rentz fac- 
tor crisis' JHenri & Sa uge 2006) have been proposed JLevinsonl 
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20071: IStern & Poutanen) 120081; iGhisellini & Tavecchiol 120081; 



Katarzvnski et alj 120081 : iBoutelier et alj 120081 : iLvutikov & Listed 
20101) . considering either significant deceleration of the inner jet 
or a multi-zone structure of the emitting region. 

To our knowledge, there have been no new reports of TeV flare 
with few- minute timescales , besid es the July 2006 outburst of PKS 
2155-304 jAharonian et alj|2007l) and the June- July 2005 events 
in Mrk 501 JAlbert et alj|2007[). However, there are mor e details 
known about the former event. lAbramowski et alj ( 120 101) showed 
that the high-activity state of PKS 2155-304 lasted for only a few 
days and consisted of many closely following 5-10 min long flares. 
This was a truly exceptional event, as during the remainder of the 
years 2005-2007 the source was quiet and remarkably stable. The 
TeV flux distribution is consistent with a superposition of two log- 
normal distributions of different spectral behaviour, indicating that 
the underlying mechanism of the flare is different from the one 
responsible for the quiescent state. Simultaneous X-ray and op- 
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deal monitoring JAharonian et alj|2009bl) indicates strong Comp- 
ton dominance and an almost cubic correlation between TeV and 
X-ray fluxes. This indicates that soft emission from the dominant 
gamma-ray producing region is swamped by radiation produced in 
a different region. A viable model of extreme TeV flares should ac- 
count for their apparently very low duty cycle and very high Comp- 
ton/synchrotron luminosity ratio. 

In previous works, the minijet emitting region has been treated 
as a compact blob propagating through the main jet. Here, we inves- 
tigate a more physically motivated model. We consider stationary 
outflows fueled by steady relativistic Petschek reconnection in one 
or more locations of the main jet. As the fluid is eventually slowed 
down by a terminal shock separating it from the so-called magnetic 
islands, two emitting regions that are relativistically boosted with 
respect to each other form, enhancing the relative importance of 
the inverse Compton (IC) process. In addition, if many such mini- 
jets f orm in closely aligned pairs with opposite propagation direc- 
tions JGiannios et al.ll201Cf) . it is plausible that these emitters will 
interact radiatively, leading to even greater Compton dominance. 
We study the effects of the presence of a weak guide field, radia- 
tive cooling and Compton drag. Given the magnetization parame- 
ter required to obtain high enough bulk Lorentz factors, a ~ 100, 
magnetic reconnection will accelerate electrons to random Lorentz 
factors 7 ~ 10 4 , which puts the Klein-Nishina limit for observed 
photon energy in the TeV band. Thus it is also important to care- 
fully model Klein-Nishina effects on TeV emissivity. 

We will introduce a physical scenario leading to the emer- 
gence of the minijets in Section [2] the geometrical and dynamical 
structure of a single minijet in Section[3]and a model for broad-band 
emission in both analytical and numerical approaches in Section|4] 
Our results are presented in Section__] followed by a discussion in 
Section[6]and a summary in Section|7] 
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Figure 1. Schematic representation of a toroidal current sheet at the bound- 
ary of two magnetic domains of opposite polarity, showing orientations of 
magnetic fields (blue), currents (red) and minijet flows (black arrows) from 
X-points to O-points. 



They are observed a s much more relativistic flows than the jet itself 
JGiannios et alj2009l) . 

Alternatively, the minijets could arise without the necessity 
for magnetic polarity reversal, if the Poynting-flux-dominated jet 
undergoes kink instabili ties (e.g., lEichlen 1 19931 : iBegelmanl 1 19981 ; 
iGiannios & SpruifeOOq) . Twisted magnetic flux tubes would form 
loops when they come into contact and magnetic field lines carried 
by colliding sections would be locally inverted. A current sheet will 
then form along the original contact plane and a pair of minijets 
will be perpendicular to the tube axes in the frame co-moving with 
the flux loop. The tube axes need to be very accurately aligned, 
otherwise the reconnection will proceed with a significant guide 
field, which, as we show later, is not a favourable situation. 



2 THE ORIGIN OF MINIJETS 

Magnetic fields dominate th e energy flux in the inner part of 
AGN jets JSikora et alj|2005l) and hence they preserve an order 
imprinted at the launching region. To launch a relativistic out- 
flow from an accreting black hole s ystem, a stro ng poloidal mag- 
netic field component is required dSpruitluOlOl) . It could be ei- 
ther produced locally within an ac cretion disc by the Parker in- 
stability (e.g., [Tout & Pringle 1992) o r advected from t he galactic 
environment (e.g., Spruit & Uzdensky 2005; Igumenshchev 2008; 
iRothstein & Lovelace! 2008; Bec kwith et alj|2009h . Both scenarios 
can naturally produce separate regions of opposite magnetic polar- 
ity. Their topolog y will be preserved in to the innermost region of 
the accetion disk beckwith et alj|2008l) . 

In the simplest configuration, two such domains are separated 
radially in the accretion disk and then consec utively pass through 
the jet-launching region dLovelace et al.ll 19971) . This results in two 
magnetic domains in the inner jet separated by a current sheet 
dSikora et al.l 120031) . As the magnetic fields of either domain are 
quickly stretched in the toroid al direction by lateral expansion of 
the jet teegelman et al.lll984t) . the current sheet can approach a 
toroidal, disk-like shape. Due to tearing mode instability, it sepa- 
rates along the azimuthal angle into a sequence of X-points and 
O-points (see Fig. [I}. Minijets form as coherent flows from the X- 
points to the O-points. If the magnetization parameter of the jet 
is large, the minijets are relativistic in the frame co-moving with 
the current-sheet. They are directed mainly in the toroidal direc- 
tion, i.e. perpendicularity to the jet propagation in the jet rest frame. 



3 THE STRUCTURE OF A SINGLE MINIJET 

To describe the properties of a mi nijet, we adopt t he scenario of 
relativistic Petschek recon nection i Petschek 1 19641) with a guide 
field studied analytic ally bv lLvubarskvl 1 120051) and numerically by 
IZenitani et alj J2009I) . In the jet co-moving frame the current sheet 
is in the :ry-plane, with the minijet outflow along a;-axis (see Fig. 
|2j- It is assumed that interaction between fast reconnected plasma 
and slow magnetic island leads to a stationary shock located at 
some x = fa. This shock separates what we define as the minijet 
region from the island region. Parameters describing the jet flow 
are denoted with subscript '1', those measured in the minijet re- 
gion with subscript '2' and those measured in the island region with 
subscript '3'. The reference frame co-moving with the jet fluid is 
denoted by Ox, the one co-moving with the minijet fluid by 02 and 
the one co-moving with the island fluid by O3. Quantities denoted 
by " are measured in O2 and those with '" in O3 . 

Jet plasma is assumed to be cold and highly magnetized 
(<Ti ~ 100). The magnetic field contains the antiparallel recon- 
necting component (within the iz-plane) of strength B\ and the 
guide field (along t he ?/-axis) of stre ngth Bi,g = aB\. Accord- 
ing to the model of |LyubarskyJ ((2005J), the reconnection outflow is 
thermal-pressure-dominated for a < l/(2y / oT) and magnetically 
dominated otherwise. Although this limiting value of a is not dy- 
namically important, it has a profound influence on the radiative 
properties of the system. In the following, we will consider mod- 
els with no guide field (case I, a = 0) and models with a weak 
guide field (case II, a — l/(2y / o r T)). The inflow into the reconnec- 
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Figure 2. The geometric setup of our minijet model. The minijet region 
{left) is separated from the island region (right) by a stationary shock (yel- 
low). Reconnected plasma flows along x axis. Spatial dimensions of the 
model are determined by the length I2 and the opening angle 1^2 of the 
minijet region. 



tion region is not relativistic, thus the relation between density and 
mag netic field i s pic 2 ~ Bf/(4-K<Ji). 

iLvubarskvl 020051) estimated the parameters of the reconnec- 
tion outflow: bulk Lorentz factor F2 ~ ^Jcn, density p2 ~ 



oTpi. The reconnected magnetic field is B2 ~ BiOif-y/ai, 
where #1 is the angle between the magnetic field lines and the 
oblique shock surface in the inflow region. However, if a guide field 
is present, it will be compressed to the value B2.G ~ 2 v / <7i\Bi,g- 
In case II this yields i?2,G ~ Bi, greatly exceeding the re- 
connected component. Pressure in the minijet region is given by 
P2 = 7 3 2 /(87r) in case I. As demonstrated in numerical simula- 
tions by IZenitani et al.l J20091) . this parameter is very sensitive to 
the strength of the guide field. We adopt the following scaling: 
P2 ~ (1 — 2(Tia 2 )_B 2 /(87r), thus in case II the minijet pressure 
is roughly half of the value in case I. The ratio of thermal to rest 
energy densities is P2/{p2<?) > y/cFi/8, thus minijet matter is rel- 
ativistically hot for <xi > 60. The magnetization of the minijet is 
02 ~ f? 2 /(2(Ti) in case I and U2 ~ 1 in case II. The minijet has an 
opening angle ifa ~ f?i/(2ai) and volume V2 ~ 1P2I2, assuming 
that its width Ay is similar to its length I2. 

The propagation of the minijet flow is affected by radiative 
processes described in Section 14.2.11 radiative losses of the elec- 
trons tend to reduce their pressure, while Compton drag by pho- 
tons emitted from the island region causes deceleration of the bulk 
flow. Thus, we introduce correction factors that describe the final 
values of evolving parameters: P2S — r\pP2 and T2,f = "T)y^2- We 
will calculate these factors using the numerical scheme described 
in Section|43l 

Initial parameters of the island region can be found by solv- 
ing the shock-jump conditions. Under the assumption of negligible 
matter rest energy density on both sides of the shock and highly rel- 
ativistic upstream flow, one can reduce the problem to a quadratic 
equation for F3, giving the following solution (see Eq. 4.11 in 
iKennel & Coronitilll984l) : 

r 3 ~ 7^7-1 1 , , N \%°h + 26o- 2 , f + 17 + 



16(a 2 , f + 1) 

+ (2ov + r 



16cr 2 f + 16t7 2 ,f + 1 



where o"2,f = B%/(16wP2,f). Then we find: 

0"2,f 



(T.i 



B, 



P3 



(cr 2 ,f + l)/3 3 - 02,f ' 

r 2 ,f 



r 3 /V 
r a , f 

r 3/ 3 : 



-B 2 



■P2. 



(1) 

(2) 
(3) 
(4) 



167T(T3 

Since Ts <C 1*2, radiation from the island region is strongly 
boosted in the minijet co-moving frame. And we expect that due 
to plasma compression both the magnetic field and average parti- 
cle energy may be significantly higher than in the minijet region. 
We assume that the length of the island region is of the order of its 
height I3 ~ 2ip2^2, hence its volume is V3 ~ Wl- 



4 EMISSION MODEL 

4.1 Particle energy distribution 

Non-thermal particle acceleration in relativistic plasmas is an 
open field of research. There are numerous studies of this 
process in relativistic shocks {e.g., Bednarz & Ostrowski 1998; 
Achterberg et al.l200ll;lLemoine & Pelle tier 2003; Nishika wa"et al.l 



20031 ; ISpitkovskvll2008l ; Isironi & Spitkovskvl 120091) . leading to a 
self-consistent (at least in weakly magnetized pair plasmas) sce- 
nario in which the first order Fermi process is initiated by the 
Weibel instability, producing power-law particle spectra N oc 7~ p 
of index p ~ 2 — 3. Particle acceleration in current sheets un- 
dergoing relativistic magnetic reconnection is in principle a more 
straightforward proc ess, since strong ele c tric fields are present . 
Numerical studies ([Larrabee et al.l 120031: Uaroschek et alj I2004J ; 



IZenitani & Hoshinoll200l 120081 ; iLvubarskv & Livertsll2008l) have 
shown that a relatively hard particle spectrum (p ~ 1) can be pro- 
duced, if the current sheet can be stabilized against relativistic drift- 
kink instabilities (RDKI) by introduction of the guide component 
of magnetic field. Particles energized to the point that they are able 
to leave the current sheet can be further accelerated by the Fermi 
process , as they bounc e between the two regions of reconnecting 
inflow dGianniosll20irj|) . 

We are not investigating details of the particle acceleration 
process here, but consider the injection of relativistic electrons of 
fixed energy distribution. In the case of PKS 2155-304, the electron 
distribution cannot be constrained by multiwavelength observations 
if we accept the argument that X-ray emission is not produced co- 
spatially with the TeV emission. However, the H.E.S.S. results im- 
ply the existence of a soft high-energy t ail that in the flaring stat e 
shows a harder-when-brighter behaviour JAbramowski et alj2010h . 

The dynamical model of relativistic reconnection and the 
shock jump conditions provide us with the values of density p and 
pressure P in both the minijet and post-shock region. We assume 
that internal energy is equally divided between protons and elec- 
trons. The average energy of an electron depends on plasma com- 
position and is highest when no electron-positron pairs are present. 
The electron number density is then given by n e ~ p/m p . The 
equation of state for both particle species is P e Cp) = 5e(p)e e ( P )/3, 
where e is the internal energy density and g is a parameter equal to 
1 for relativistic particles and 2 for non-relativistic particles. From 
equipartition we have e e = e p , hence the pressure of electrons 
(protons) is 



e(p) 



9e(p) 

le + g. 



while their average random Lorentz factor is 



<7c( P )) ~1 + 



e c ( P ) 



n e(p) TO e(p) ( 



1 + 



3 



pc 2 g e + g P m e ( P ) 



(6) 



(7) 



It is clear that electrons are highly and protons at least mildly rela- 
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tivistic, since P/(pc 2 ) > 1. This points to a strongly peaked elec- 
tron energy distribution. 

For simplicity, we will adopt a relativistic Maxwellian electron 
distribution 



nc(7) = nc^e" 7/7 ° 
27o 



(8) 



where 70 = (7c) /3. In models intended to fit the observed spectra, 
we add a power-law high-energy tail 



Tie (7) = n c 



27 J 

7~ p 



P+2NP+2 



(**) 



7 < (p + 2)70 
-7/7cut otherwise 



where p > 2 is the electron index, 



2- 



7o = <7c) 



[' 



+ 1 



, + W+ (£+?f\ t 



6 - [p 3 + 9p 2 + 30p + 38 + {j ^f\ 



-(p+2) 



(9) 



(10) 



and 7cut 2> (p + 2)70 is the Lorentz factor of the exponen- 
tial cut-off. For this choice of normalization factors and 70, the 
two distributions and their first derivatives are joined continuously. 
A different combination of relativistic Maxwellian and power- 
law components of the electron distribution has been proposed by 
iGiannios & Spifkovskvl ( 120091) . 



4.2 Radiative processes 

In this Section, we calculate emission from the combined mini- 
jet and island regions. To this end, we consider synchrotron radia- 
tion, IC scattering (with Klein-Nishina cross-section) of local syn- 
chrotron photons and of radiation from the other emitting regions; 
radiative losses from all processes; Compton drag from anisotropic 
photon fields in the minijet co-moving frame and 77 pair produc- 
tion absorption. In Table [TJ we compare parameter values in both 

cases I and II for reference values o\ — 100, B\ = 10 G, 

I2 — 10 14 cm and Tj c t = 10 taken from lGiannios et al.l J2009l) . 



4.2.1 The minijet region 
Synchrotron emissivity is given by (see Eq, IA2b 
■" <JtC / 2 \ ,,,.. 

JSYN,2 = -5— n e , 2 MB,2 <7e,2> • (H) 

The photon energy of the synchrotron peak is at /ws'yn,2 ~2x 
10~ 8 (-B2/I G) (7^,2) e V. Energy density of the synchrotron 
radiation is given by (see Eq, IC2b 



""SYN.2 — 



JSYN,2 



dV 



r 2 c J v r 2 (l + (3 2 cos6) ' 



(12) 



where 9 is the angle to the emitting element measured from the 
flow velocity direction and r is the distance to the emitting vol- 
ume element. As the energy density of the synchrotron radiation 
is dominated by photons coming from the backwards direction, 
it should depend strongly on the position in the minijet region. 
We can estimate it, by taking into account only the solid angle 
of 8 > (it — 1/T), as Usyn,2 ~ 2-i/j2JsYN,2^2/c. Using this, 
we calculate synchrotron-self-Compton (SSC) emissivity Js'sc.a ~ 
Jsyn 2( w syn 2/ u b,2). The characteristic frequency of the SSC 
component is /i^ssc,2 ~ (7^,2) h,Ug Y N,2- Klein-Nishina effects 
become important if bggc 2 = 4 (7e,2) hvg Y N,2/( rn e c ' 2 ) ^ !• 



From TableQ] we find that SSC emissivity is ~ 3 orders of magni- 
tude weaker than synchrotron emissivity and is safely in the Thom- 
son regime. 

The energy density of radiation produced behind the stationary 
shock, in the island region, can be found from (see Eq. lClb : 



"3^2 = 






(l + /3 2 cos6>) 2 dV 



(13) 



l V3 (l + p 3 cos9)3 r 2 ■ 

This radiation is observed mostly at small 9, thus we can set 
cos 9 ~ 1 to obtain an estimate (at x — Z2/2): 



W3-»2 ~ 1671-^3 i 



33 



1(1+ /3s 



using which we find external Compton (EC) emissivity 
// 

n2 w 3->2 



JeC,2 3, „ 

- -(1 + ^obs 

JSYN.2 ^ 



MB, 2 



(14) 



(15) 



where /i" bs = (cos Sobs + /32)/(l + P2 cos0 o bs) is the cosine of 
the observer inclination in O2 and # ODS is the observer inclination 
in d (see Section |4.2,3b , For the IC process only synchrotron pho- 
tons from the island region are relevant; their characteristic energy 
in 02 is hv'^2 ~ 2^/^3(1 + /?3)]^syn,3- The characteristic 
energy of EC photons is ^0,2 ~ 2 (! + Mobs) (7*2) ^2 and the 
Klein-Nishina parameter &ec,2 ~ h^EC^/dle,?} m e c 2 ). We find 
that the EC component is stronger than the synchrotron component 
by ~ 4 orders of magniture in case I and 1 order of magnitude in 
case II. However, characteristic photon energies lie ~ 2 orders of 
magnitude within the Klein-Nishina regime and EC peaks should 
be suppressed. This means that in case II, synchrotron is likely to 
dominate. 

From the relative importance of the spectral components we 
find that the radiative cooling of electrons is dominated by comp- 
tonization of the island photons in case I and by synchrotron emis- 
sion in case II. Using Eq. l lAllb . we estimate the efficiency of elec- 
tron cooling for plasma crossing the minijet length: 



Ccool,2 — — 



d7 e h 



^CF T (7e,2) 



" 7 

u 2 L 2 



(16) 



dx (7e,2) 3m e c 2 

where w 2 stands for the largest of Ug 2 an d %xt,2- As shown 
in Table [TJ we obtain large cooling efficiency in case I, when 
u ext,2 ^S> «b ,2- However, as we are in the Klein-Nishina regime, 
this result is overestimated. In case II the cooling is not very effi- 
cient. 

When electron cooling in the minijet region is dominated by 
IC losses off post-shock radiation, highly anisotropic in the minijet 
co-moving frame, and the minijet plasma is pressure-dominated, 
Compton drag may result in significant bulk braking of the minijet 
flow. The gradient of bulk Lorentz factor is given by (see Eq. 12 in 
ISikoraet all 19961) : 



dr 2 

dx 



drag 



M c ,2 + M Pi 2 + MB, 2 



(17) 



where u, 



o(p),2 



is the total co-moving energy density of electrons 
(protons), ub, 2 = Sf/(8 7r ) I s me energy density of magnetic field 
and p'e^il) is me average co-moving radiative force per electron 
of random Lorentz factor 7. The full Klein-Nishina formula for 
p e (7) is provided in Appendix iBl (Eq , lB3l . In the Thomson regime 
it becomes p" 2 (7) ~ (2/3)ctt7 2 M3_ +2 - We estimate the efficiency 
of the Compton drag as: 



^,drag,2 



dTa h_ _ 2 (7^,2) (n e ,20- T l2)u'^2 

' dx r 2 ~ 3 T 2 p 2 C 2 + 3P 2 + M B ,2 



(18) 
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We find that radiative drag should be more efficient in case I than 
in case II, however our estimate is again subject to Klein-Nishina 
suppression. We note that the SSC process can counter the radiation 
drag, as the local synchrotron radiation field is not isotropic in 02 
and is actually stronger when arriving backwards with respect to 
fluid motion (see Eg, 112b. However, since external radiation domi- 
nates local synchrotron radiation, this effect is not important in our 
problem. We left the detailed solution of this problem to our nu- 
merical scheme. 

As a consequence of the deceleration of the minijet flow, 
plasma will be compressed in its co-moving frame. We adopt here 
adiabatic scaling of plasma parameters, noting that various plasma 
effects may complicate the relations among them. For a small de- 
crease of Lorentz factor by 8Y2, number density of particles should 
increase by 8n,2 ~ (5r2/F2)n2. As the magnetic field is perpen- 
dicular, it increases like matter density 8B2 ~ [8Y2/Y2)B2- Pres- 
sure dominates the enthalpy and thus it should increase quadrat- 
ically, 8P2 ~ 2{8T2/Y2)P2- Such simultaneous changes in 712 
and P2 imply an increase of energy of each electron by <57 c ,2 ~ 
(<5r2/T2)7c,2- Note that the electrons do not gain energy in the 
Oi frame, rather a part of their bulk kinetic energy is converted to 
random kinetic energy. 

High-energy photons will be absorbed in photon-photon 
pair production mainly on the synchrotron radiation. The 
cross-section given in Eq. dA13b peaks at <j 77 ~ ctt/5 



for photon energy hv. 



77 ,peak 



3.5(r 



l /{hv t 



argct } 



In 



0i, the minijet synchrotron photons provide targets of en- 



ergy hvts. 



2V 2 hv 



SYN.2 



(as they come mainly from 



the backwards direction) and energy density Wtarget,2 ~ 
4r|-Ug Y N 2> while the island synchrotron photons have higher en- 
ergy /i^targct,3 ~ /i^syn, 3/^3(1 + ft )] and energy density 
Mtarget,3 ~ u 3-»-2 /(4r| ). Using Eq. dA12t . we estimate the mean 
free path for high-energy photons 

x 15 /^target, 2(3) 



Vy,2(3) 



(19) 



K 77 ,2(3) °T ""target, 2(3) 

In Table Q] we find, that A 77i 2(3) 3> h, thus we expect little ab- 
sorption from both sources of synchrotron photons. However, when 
increasing the source characteristic size or parameters governing 
the synchrotron emissivity, we expect this effect may put impor- 
tant constraints on VHE luminosity. Detailed calculation will be 
performed in our numerical scheme. 

4.2.2 The island region 

Increased magnetic field and pressure will produce higher syn- 
chrotron emissivity than in the minijet region. Due to the assumed 
geometry of the region, synchrotron energy density can be esti- 
mated as Usyn,3 ~ 7r isYN,3^3/(2r'3c). Energy density of radia- 
tion from the minijet region is given by 



^2^3 



d 


■// 

.h_ 
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(1 


+ P3 


cos6>) 2 


dV 


r| 


(1 


+ /3 2 


cos#) 3 


r 2 



(20) 



In contrast to the situation in the minijet region, we have mainly 
6 ~ 7r. Within the solid angle 6 > n — I/T2, we can set cos 6 ~ 
— 02, obtaining an estimate 
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(21) 



The characteristic energy of photons from the minijet region is 
hv'iUz ~ 2r 2 /[r 3 (l + /3 3 )] hVsYK,2- EC emissivity differs by 
the sign of the cosine of scattering angle: 



Table 1. Parameters of plasma composition, electron distribution and radia- 
tive output from the minijet (subscript 2) and island (subscript 3) regions 
calculated for CTl = 100, B 1 = 10 G, l 2 = 10 14 cm, 0i = 0.5, r\ P = 1, 
t)y = 1, #obs = 0.1 and Tjct = 10. Model with no guide field (case I) is 
compared to a model with a weak guide field (case II). 
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(22) 



where ^" bs = (^ obs + ft)/(l + ftMobs). The characteristic en- 
ergy of EC photons is l4'c,3 ~ 2(1 - /i" bs ) (7^3) v' 2 'U%. Other 
parameters are found in the same way, as for the minijet region. We 
find that the strongest spectral components are synchrotron radia- 
tion and SSC, while the EC component is negligible. Synchrotron 
emissivity is similar in both cases, however SSC is much stronger 
in case I. The energies of SSC photons lie ~ 2 orders of magnitude 
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into Klein-Nishina regime, similarity to EC photons in the minijet 
region. 

The radiative cooling efficiency is calculated from 



t _ d7 e h Aa T (7.5.3) in, 

Cc ° o1 ' 3 - dx ( 7e , 3 ) ~ 3m e c 2 T 3 " 3 h 



(23) 



In case I cooling would be dominated by SSC (£cooi,3 ~ 1.7), how- 
ever Klein-Nishina suppression is likely to lead to the domination 
of synchrotron cooling (Ccoo1,syn,3 ~ 0.26). In case II synchrotron 
cooling is significant, so the electrons can radiate the bulk of their 
energy within the distance I3. 



4.2.3 Obserx'ed luminosity 

Let us place the observer in the xz plane (see Fig. f2j at an- 
gle Sobs to the minijet velocity direction as measured in jet rest- 
frame Oi, so that fcobs = [cos Sobs, 0, sin Sobs]- The Doppler fac- 
tor between the minijet co-moving frame 02 and 0i is 2?2->i = 
[TMl — P2 cos Sobs)] 1 , between the island co-moving frame 03 
to 0i is ©3_»i = [r3(l — /?3 cosSobs)] 1 and between 0i and the 
laboratory frame is T>i^o — rj et (l + /3j e t sinS b s )- The luminos- 
ity seen by the observer from the minijet (island) region is 



£obs,2(3) = 47Tl> 1 _ ) . I> 2 (3 ) _nj2 0'ai")V2 



(3) 



(24) 



In order for the island region to compete with the minijet region in 
observed luminosity, the emissivity ratio should be 






>w 



(25) 



ft VUl ^3 \l3/ 4^2 

As can be seen in Table Q] this is the case for the ratio of syn- 
chrotron emissivities, since n e , ub and (7,.) are significantly higher 
in the island region. We find that the low-energy component from 
the island region should dominate in case I, while both components 
should be comparable in case II. However, high-energy emission 
is clearly dominated by Comptonization in the minijet region of 
synchrotron photons from the island region. Not only are the dif- 
ferent spectral components produced by different mechanisms, but 
they also actually come from different (albeit adjacent) parts of the 
complex reconnecting system. 



4.2.4 Variability timescale 

The minijet length is constrained by the observed variability 
timescale. The 7-ray emission is dominated by the contribution 
from the minijet region, which is static in the jet co-moving frame 
0i . The variability timescale resulting from light-travel effects cal- 
culated in this frame should be divided by the Doppler factor of this 
frame with respect to observer, X>i^o ~ Tjot. The time delay be- 
tween signals emitted by the same plasma portion at x — and 
x = h is At' x = (^2/c)(/3 2 ~ 1 — cos Sobs), which for S bs ~ T^ 1 
is of the order of At' x ~ fa/(r|c). This would be larger if mini- 
jet plasma deceleration were taken into account. Time delays be- 
tween signals emitted simultaneously in 0i across the minijet re- 
gion width (Ay ~ h) are At' y = (h/c) sinSob a ; for S bs ~ T^ 1 
it is At'y ~ h/(^2c). Unless S bs < T^ 2 , which is an unlikely 
situation, At' y 3> At' x and the minijet length can be estimated 
as I2 < rj ct r2Ct var , which is consistent with the estimate for 



a sphencal blob prop a 
JGiannios et al]2009l) . 



4.3 Numerical scheme 

We use exact formulae summarized in AppendixlAlto calculate the 
spectra of radiation emitted from the minijet and island reg ions. 
We use elements of the BLAZAR code JModerski et alj|2003l) wit h 
formulae valid in the Klein-Nishina regime dModerski et alj2005h . 
The minijet region is divided along the x axis into N sec- 
tors of width (Ax) 1 , median position x l and volume AV 1 = 
2ip2h (Ai)' x 1 fixed in 0i. For each sector, the synchrotron emis- 
sivity spectrum (j'syn 2YW) is calculated from Eq. dAlt . syn- 
chrotron energy density (wg Y N 2Y( V ") from Eq. | |12I >, SSC emis- 
sivity (j'iscfiYiy") from Eq. dA7t . energy density of radiation 
from the island region 



(«3->2)V) 



T-l •III/ lll\ 

T2 n \v ) 



r 2 
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1 + fh cos S dV 
V3 (l + ftcosS) 2 ^ 2 " 



(26) 



(see Eq. I13t and EC emissivity (Jec,2Y{ v ") from Eq. !A3t , En- 
ergy densities of synchrotron radiation and the radiation from the 
island region are taken into account when calculating electron cool- 
ing rate, Compton drag and photon absorption. 

Compton drag efficiency is calculated from Eqs. d!7| > and d!8l >: 



Cdrag,2 — 
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dT| 



(27) 



drag 



The evolution of minijet fluid parameters due to the plasma decel- 
eration and compression is 



r 2 -* 1 - C 



^drag,2 I ^2 > 



(28) 
(29) 
(30) 
(31) 



The electron distribution is initialized to n ei 2 in the first sector, 
then both effects of adiabatic compression and radiative cooling are 
taken into account in their evolution: 



bi^ [1 + a ag , 2 ) bi , 

'A -*■ (l + Cdrag,2j ™2 ' 
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dx 



where we employ a kinetic equation 



dn e , 2 (7) 
dx 



9 ( . . d 1 
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(32) 
(33) 

(34) 



and j^|rad is the cooling rate given by Eq. dAllt . 

The evolution is tempered by including the effect of mixing 
the evolved fluid with freshly reconnected plasma proportionally to 
their volumes: 



XI 



(AxY/2] xXJ + (Ax)' x X 2 
x i + (AxY/2 



(35) 



where X2 stands for r 2 , B2, JI2, P2 or n e ,2(7). 

Essentially, at every step we mix x' — (AxY/2 parts of old 
electrons evolved via radiative losses with (AxY parts of new elec- 
trons injected through the minijet boundary. The nu merical method 
for so lving the kinetic equation is explained in iModerski et al.l 
J2003h . 

From the final electron distribution n^ we find (7e,2) , 
and then the pressure 
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DiV + l 



<7e,2) 



JV+1 



"5e,2 + g v ,2 



(36) 



<7e,2) '' J .9e,2 + 3p,2 

Our approach to properly accounting for the effects of electron 
cooling and Compton drag is to find such values of corrections rjp 
and r/r, for which P 2 Ar+1 = P 2 ,f and T^ +1 = I\ t . 

In a similar way we calculate emissivities in the island region, 
which is also divided into several sectors, however the effects of 
Compton drag and plasma mixing are not considered in this case. 
We then calculate the luminosities from each sector and transform 
them to the laboratory frame, using Eq. ( 124b . 

Pair production opacity is calculated in 0\ by following a 
ray emitted at the center of each sector in several steps evenly 
spaced over distance 2^2- The minijet region is divided into a 2- 
dimensional array (j,k) £ (1, ..., N) x (1, ..., N) of sectors cen- 
tered at x j = (j - l/2)l 2 /N and y k = [(k - l/2)/N - l/2]/ 2 
of volume AV 2 3 ' = 2^2 ^1 aj "' /A'" 2 , while the island region is di- 
vided into a 1 -dimensional array (k) £ (1, ..., TV) of sectors cen- 
tered at x = h + h/2 and y k = [(k - l/2)/N - 1/2]Z 2 of 
volume AV3 = 21P2I2I3/N. For each volume element, the mean 
distance r mea n from the photon, mean scattering angle and mean 
Doppler factor of the emitting fluid are calculated and used to de- 
termine energy density of incident radiation at given position. A 
problem arises when r mea n is small or comparable to the size 7? 
of the volume element. For a spherical volume element contain- 
ing static, isotropic and optically thin emitting fluid, the radiation 
density at distance r mcan from the center is given by 



u(r n 



3X 

tvR 2 c 
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2q 
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(37) 



where q = r mcs . n /R and L = 4irjV(R) is the total luminosity. 
This formula can be well approximated for both small and large 
values of q with the following: 

L 



U\Tmc&n) 



(38) 



We thus replace the mean distance with r = r mcan + R /3 and 
use formulae valid in the limit of large q. We adopt R ~ h/(2N) 
for elements of both the minijet and island regions. 



5 RESULTS 

We begin by running our numerical scheme with the same param- 
eters as used in Table [TJ We denote this model as Model A. Since 
effects like opacity, Compton drag and radiative cooling are weak, 
a direct comparison of emissivities, energy densities and luminos- 
ity components between analytical estimates and numerical results 
can be performed. 

In Fig.[3]we present spectral energy distributions (SEDs) in the 
observer's frame, electron energy distributions and pair production 
opacity as a function of emitted frequency. There is good agree- 
ment between estimated and numerically calculated positions and 
values of peak radiation energy densities, local emissivities and ob- 
served luminosities. Synchrotron emission from the island region 
has been overestimated, because effective cooling of electrons re- 
duces the average emissivity. The EC component from the mini- 
jet region has been underestimated, because emissivity from the 
final minijet sector rapidly increases in the vicinity of the station- 
ary shock front. Cooling efficiency is low for the minijet electrons 
and high for the island electrons. The pair production opacity for 
TeV photons is at a comparable level of r ~ 0.01 in both cases and 
is very uniform when measured in different geometrical sectors. 
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Figure 3. Results obtained for Model A - parameter values are cr\ = 100, 
I2 = 10 14 cm and B± = 10 G. The observer is located at 9 ^ s = 0.1. 
Red lines denote model with no guide field (case I), while blue lines denote 
model with weak guide field (case II). Analytical estimates from TablefT|are 
marked with crosses (for the minijet region) and circles (for the island re- 
gion). Top panel: SEDs of the minijet (thin solid lines) and the island (thin 
dashed lines) regions, as well as their sums (thick solid lines), in the labo- 
ratory frame. Middle panel: electron number density energy distributions 
in local co-moving frames. Evolution of the electron distribution is illus- 
trated by plotting results for several sectors out of 10 equally spaced along 
x-axis. For the minijet region (solid lines) we show sectors 2 (thick line), 4, 
6, 8 and 10 (lines are barely distinguishable as a result of inefficient cool- 
ing). For the island region (dashed lines), we show sectors 1 (thick lines), 
2, 3 and 4. Bottom panel: pair-production opacity for high-energy photons 
emitted from the minijet region contributed by soft photons from the minijet 
(solid lines) and island (dashed lines) regions. We plot opacities for photons 
emitted from sectors 2, 4, 6, 8 and 10 of the minijet region. The results for 
different sectors are barely distinguishable, thus we use the same line types. 



It is dominated by the contribution from the island synchrotron 
photons, since absorption by minijet synchrotron photons targets 
higher-energy emission. This is in qualitative agreement with tar- 
get energies and mean-free paths listed in TableQ] 

In Fig. [4] we show radiation energy densities and emissivities 
in local frames. There is strong evolution of the synchrotron emis- 
sivity from the island region (solid lines in bottom right panel), 
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Figure 4. Local energy densities and emissivities for model A - same as in Fig. [5] Analytical estimates from Table[T]are marked with symbols explained below. 
Left panels: the minijet region - sectors 2 (thick lines), 4, 6, 8 and 10 out of 10. Right panels: the island region - sectors 1 (thick lines), 3, 5, 7 and 9 out of 
10. Top panels: energy density of the magnetic field (solid lines and crosses), synchrotron radiation (dotted lines and circles) and external radiation (dashed 
lines and squares). Bottom panels: synchrotron (solid lines and crosses), SSC (dotted lines and circles) and EC (dashed lines and squares) emissivities, all 
calculated in their respective local frames. Line colours have the same meaning as in Fig. [3] 



reflected in the local synchrotron radiation energy density in the 
island region (dotted lines in top right panel) and the external radi- 
ation energy density in the minijet region (dashed lines in top left 
panel). This, however, does not affect the dominant high-energy 
emissivity components: EC in the minijet region (dashed lines in 
bottom left panel) and SSC in the island region (dotted lines in bot- 
tom right panel), which agree with the estimates. 

In Fig. [5] we show the effect of varying the observer's incli- 
nation angle o bs on observed SED and opacity. For the values of 
this angle of the order of a few times I/T2 ~ 0.1, the Doppler fac- 
tor of the minijet region X>2^i varies much more than the Doppler 
factor of the island region 2?3_>i. The flux produced in the minijet 
region, which dominates in soft X-ray and GeV-TeV bands, de- 
creases strongly with increasing 8 b B - On the contrary, the hard X- 
ray flux originating in the island region slightly increases, because 
the Doppler factor of the jet £>i->o increases faster than £>3^i de- 
creases. For # ba ^ 0.5, emission from the island region dominates 
the total output. Pair-production opacity shows variation of thresh- 
old energy with increasing 8 b s '■ a moderate decrease for the contri- 
bution from the minijet photons and a small increase for the contri- 
bution from the island photons. However, the contribution from the 
island photons always dominates and the net effect for TeV photons 
is a small decrease of opacity with increasing (9 bs- 



5.1 Maximum TeV luminosity 

iGiannios et al.l J2009T) chose the parameter values for Model A in 



observed fast TeV flares. The emitting region size I2 



10 1 



was estimated from energetic considerations corresponding to a jet 



magnetic field value of B\ ~ 12 G, while the variability timescale 
constraint was h < 9 x 10 14 cm for i va r = 5 min. We find that 7- 
ray luminosities for Model A are of the order of 10 43 erg s _1 , ~ 4 
orders of magnitude lower than required. This discrepancy likely 
comes from the differ ent geometrical shape of the emitting region. 
IGiannios et all J2009I) used a highly idealized blob model for the 
emitting region, but a more realistic description of relativistic mag- 
netic reconnection requires extremely small opening angles of the 
emitting regions, which limits the emitting volume. To increase the 
luminosity, we could increase either the energy density regulated 
by B\ and/or the minijet length. We will therefore consider mod- 
els with a minijet length approaching the maximum allowed by the 
variability timescale and tune the value of the jet magnetic field to 
produce the desired value of the luminosity. 

It is then interesting to determine the maximum TeV lumi- 
nosity that can be produced by minijets and to identify the physi- 
cal effect that limits it. An obvious limitation on 7-ray luminosity 
comes from the intrinsic opacity to high-energy radiation. Another 
constraint comes from radiative efficiency. If the electrons are not 
heated everywhere accross the source, the effective emitting vol- 
ume will depend on the average cooling distance of high-energy 
electrons. In our model, efficient Compton drag could also be a fac- 
tor limiting high-energy emission, as it decreases the boosting fac- 
tor for part of the emitting volume. All these effects become more 
prominent with increasing energy density of soft radiation and our 
numerical scheme is designed to deal with such a problem. 

The relativistic Petschek reconnection model from iLvubarskvl 
120050 relates particle density and pressure to the magnetic field 
strength, thus we can increase radiative energy densities and emis- 
sivities by changing one parameter, without altering the geometric 
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Figure 5. The effect of changing the observer's inclination angle 9 b s . as 
measured with respect to the flow velocity (x-axis) in the jet co-moving 
frame Ox, on the observed SED (top panel; summed emission from the 
minijet and island regions) and pair production opacity (bottom panel; for 
the minijet photons emitted from sector 10; by all soft photons) for Model 
A. The values of 6* obs are (solid lines), 0. 1 (default value for other mod- 
els), ..., 0.6 (dashed lines). Line colours have the same meaning as in Fig. 

m 



structure. In Fig. [6] we show a series of models with increasing jet 
magnetic field. We find that on average, an increase of magnetic 
field by a factor 2 increases both synchrotron and IC luminosities 
by a factor ~ 11 and opacity at v — 10 Hz by a factor ~ 5. 



The opacity approaches unity when high-energy luminosity is at 
the level of ifj 47-48 erg s -1 . Cases I and II do not differ signifi- 
cantly in terms of luminosity (with the exception of the soft X-ray 
synchrotron peak from the minijet region) and opacity, however, 
Compton drag is more efficient in case I, as already predicted in 
Section |4.2. 11 For the models with opacity close to unity, the bulk 
Lorentz factor in the minijet region drops by ~ 60% in case I and 
only by ~ 5% in case II. Thus, the Compton drag effect can be 
an important factor limiting high-energy luminosity in Compton- 
dominated anisotropic sources. 

The conclusion from the previous paragraph is that for a mini- 
jet size of 10 14 cm, jet magnetization ax — 100 and bulk Lorentz 
factor Tjct = 10, the 7-ray luminosity is limited by opacity at 
~ 10 48 erg s" 1 . Can we relax this constraint by considering the 
Comptonization of external radiation? Note that our scenario sug- 
gests a chain of minijets forming a ring-like structure. Each mini- 
jet should have two neighbours, one sharing the X-point and one 
sharing the O-point. Both would be directed away from the exter- 
nal observer and thus relativistically hidden. The average inclina- 
tion angle between two neighbours is inversely proportional to the 
number -/V r i ng of individual minijets in the ring. If N r i ng > 7TT2, 
radiation from the O-point-sharing neighbour is strongly boosted 
in the comoving frame of the observed minijet region (02). It is 
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Figure 6. The effect of increasing the minijet energy density, related by 
the reconnection model to the jet magnetic field strength Bx, on the ob- 
served SED (top panel; summed emission from the minijet and island re- 
gions), pair production opacity (middle panel; all contributions summed 
for the minijet photons emitted at x = O.9/2) and minijet dynamics (bot- 
tom panel; evolution of the Lorentz factor along the minijet region). The 
models have been calculated for minijet length I2 = 9 X 10 14 cm. The 
values of B\ (in Gauss) are: in case I (red lines) - 4, 8, 16 (dashed lines) 
and 32 (solid line); in case II (blue lines) - 4, 8, 16, 32, 64 (dashed lines) 
and 128 (solid line). 



almost completely anisotropic and can dominate the Compton drag 
effect. 

In Fig. [7] we show a series of models taking into account the 
radiative interaction with an opposite minijet system. It is assumed 
that the opposite minijet has identical parameters and is also af- 
fected by interaction with observed minijet. We find that the major 
difference that this additional radiation makes is to increase the ef- 
ficiency of Compton drag. In the case of a single minijet we only 
calculated the drag effect from synchrotron photons originating in 
the island region. They had relatively high energy and the bulk of 
them were scattered in the Klein-Nishina regime. Radiation from 
an oppositely-directed minijet is highly anisotropic, regardless of 
the region from which it was emitted. Photons from the minijet re- 
gion have lower energies and are more effectively scattered, even 
though they are especially strongly boosted in O2 frame. 
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Figure 7. Same as in Fig. [6] but taking into account radiation from an op- 
posite minijet system. The values of B\ (in Gauss) in both cases are: 4, 8 
{dashed lines) and 16 (solid line). 



Strong Compton drag limits the maximum magnetic field 
strength for which our model is valid, i.e., for which the deceler- 
ation is not catastrophic. This limits the resulting X-ray luminosity 
(which is of synchrotron origin), especially in case II. Because of 
higher synchrotron emission from the minijet region, models with a 
guide field are more sensitive to the existence of an opposite minijet 
system. Comparing results for models calculated with B\ = 4 G, 
the differences between isolated and mirrored minijets are small in 
case I, while in case II we observe a stronger IC component, higher 
opacity and noticeable drag when radiation from an oppositely- 
directed minijet is present. 

For the brightest models that we have obtained, the 7-ray lu- 
minosity is of the same order of magnitude as in the case of an 
isolated minijet, ~ 10 48 erg s _1 . For models with a guide field 
(case II), this can be achieved with much lower synchrotron emis- 
sion, i.e., larger Compton dominance. This is an important advan- 
tage in light of observational properties of the best-studied fast TeV 
outburst. 



Figure 8. SEDs of the minijet models matching the TeV spectrum of flaring 
state of PKS 2155-304, compared with July 2006 simultaneous observations 
by H.E.S.S. and Chandra {thick black lines). Red lines show models with 
no guide field (case I), blue lines models with significant guide field (case 
II), solid lines models of isolated minijets and dashed lines models with 
radiation from an opposite minijet system ('OPP'). 



5.2 Application to TeV blazars 

We apply our model to a TeV flare of PKS 2155-304 ob- 
served by H. E.S.S. in July 2006, w hich reached luminosities 
10 ergs -1 JAharonian et al] 120071) . We are additionally con- 
strained b y the simultaneous Cha ndra observations in the soft X- 
ray band JAharonian et al] |2009b). Since the variability amplitude 
is much lower in the X-ray band and there is evidence for an almost 
cubic relation between TeV and X-ray fluxes, the X-ray luminosity 
of the minijet must stay below the Chandra result. 

We are attempting to match the observed TeV spectrum of 
photon index ~ 3.5 (above ~ 0.4 TeV), thus a power-law tail 
is required in the electron distribution of the minijet region. Be- 
cause this part of the IC spectrum is produced in the Klein-Nishina 
regime, there is no straightforward relation between the electron 
spectral index and the photon index. We find that the value p = 3.2 
(iV 7 oc 7~ p ) produces an adequate slope of the TeV spectrum. 

In Fig.(8] we show four models with TeV luminosity matching 
the H.E.S.S. observation of PKS 2155-304: with ('II') or without 
(T) a guide field and with ('+OPP') or without opposite mini- 
jet radiation. All models have been calculated for minijet region 
size I2 = 9 x 10 14 cm. The jet magnetic field values Bi for 
this models are 8.5 G (T), 35 G ('II'), 7 G ('I+OPP') and 8 G 
('II+OPP'). High-energy components are very similar for all mod- 
els, while low-energy components are widely diverse with peak lu- 
minosities spanning 4 orders of magnitude. More soft radiation is 
produced in models with a guide field and without opposite minijet 
radiation. The low-energy component is broader in the presence of 
a guide field, extending up to ~ 100 MeV instead of ~ 1 MeV 
without a guide field. The guide field affects also the spectral shape 
of the low-energy component in the hard X-ray band: it is sharply 
peaked in case I, while in case II it shows a soft power-law plateau. 
The constraints imposed by Chandra observations exclude model 
'IT, which severely overproduces X-ray flux. Model 'II+OPP' is 
marginally consistent with the data at the high-energy end of the 
observational band. 

As already mentioned in the previous subsection, a single cur- 
rent sheet may produce more than one minijet at the same time. 
Since there is no preference for any direction in the plane per- 
pendicular to the jet (which is true for a ring-like structure), in- 
dividual minijets would arise with different orientations and most 
of them would be strongly misaligned with respect to the line of 
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Figure 9. SEDs of a system of iV r i n g = 30 minijets evenly spaced around 
the jet axis, so that only one is directed close to the line-of-sight. Individual 
minijet spectra are shown with dashed lines, summed spectrum with a solid 
line. Red lines: minijets calculated with model including Comptonization of 
radiation from the opposite minijet for case I. Blue lines: minijets calculated 
for case II (only first 3 minijets and the sum off all 30 are shown for clarity). 



sight. As we showed in Fig. [5] radiation from the minijet is strongly 
anisotropic, but relativistic boosting is much stronger in the minijet 
region, which dominates 7-ray emission, than in the island region, 
which produces strong hard X-ray emission. Emission from mis- 
aligned minijets should be taken into account in a discussion of 
X-ray constraints. They can also contribute to the 7-ray luminosity 
when -/V r i n g ^> 7rr 2 . 

In Fig.|9]we show observed luminosities of a series of -/V r i ng = 
30 minijets evenly spaced around the jet axis, forming a ring shown 
schematically in Fig. [JJ Each minijet has exactly the same pa- 
rameters and they differ solely by the observer's orientation. Only 
one minijet is closely aligned with the observer (#obs = 0.1; the 
highest line) and every consecutive one is rotated by an angle 
27T /Wring = 0.21 in the 0\ frame. In case I (red lines), we find that 
even the second-best-aligned minijet (second highest line) is a neg- 
ligible 7-ray emitter. However, in the hard X-ray band the ensemble 
of all misaligned minijets will contribute significantly. The summed 
spectrum of all individual minijets will be only marginally con- 
sistent with Chandra limits. Therefore, the number of misaligned 
minijets per each aligned one should be less than 30 to explain 
the flare in PKS 2155-304. In case II (blue lines), however, emis- 
sion from misaligned minijets is negligible across the whole en- 
ergy range. This is because the Lorentz factor in the island region 
is much higher for this model (TI+OPP'; 1^ = 2.13) than in model 
T+OPP' (r 3 = 1.06). 

5.3 Lowering jet magnetization 

Jet magnetization of (Ti ~ 100 is required for reproducing fast TeV 
flares in the minijet model. It enables high minijet Lorentz factor r 2 
and large average electron Lorentz factor (7e,2), both of which are 
proportional to y/ai. The latter requirement can be eased, noting 
that some sort of stochastic particle acceleration process is neces- 
sary in order to reproduce power-law spectral tails to fit observa- 
tions. Lower <j\ implies lower r 2 and thus a more compact minijet 
region for the same observed variability timescale. Reduced emit- 
ting volume would have a strong impact on the luminosity of both 
spectral components. We study this effect in the class of models 
including opposite minijet radiation, trying to match the TeV lumi- 
nosity of PKS 2155-304 or, if impossible, calculating a model of 
maximum luminosity. 

In Fig. [10] we compare the SEDs for three values of o\ . The 
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Figure 10. SEDs of minijets calculated for u\ = 100 (solid lines; same 
as the dashed lines in Fig. [8}, 50 (dashed lines) and 25 (dotted lines). The 
models have been calculated for cases of no (red) or weak (blue) guide field, 
including radiation from the opposite minijet. 



value <ti = 100 (solid lines) has been used in lGiannios et alJa2009h 
and in previous paragraphs. o\ = 50 (dashed lines) corresponds to 
h = 6.4 x 10 14 cm and F 2 = 7.1, while a\ — 25 (dotted lines) 
to l 2 = 4.5 x 10 14 cm and T 2 = 5. We were able to fit H.E.S.S. 
data for PKS 2155-304 for ai = 50, but not for <ri = 25, where 
opacity limits TeV luminosity below the observed level. Keeping 
Tjot = 10, the last case corresponds to effective Lorentz factor of 
the minijet plasma r 2 Fi c t ~ 50, the minimum value derived by 
iBegelman et all ( 120081) . Thus, our model confirms that prediction, 
even though it has been derived within a single-zone framework. 



6 DISCUSSION 

Our calculations show that it is much easier to obtain a high Comp- 
ton dominance for minijet models based on relativistic magnetic 
reconnection with no guide field (case I). Inspection of Table[JJre- 
veals that this is related to two factors. First is a significantly lower 
magnetization of the minijet region plasma <r 2 , which regulates the 
ratio of magnetic to electron pressure. For roughly the same ther- 
mal energy carried by particles in both cases, the magnetic energy 
density is 2 orders of magnitude lower in case I, and so is the syn- 
chrotron emissivity. The second reason is the much stronger com- 
pression of plasma crossing the stationary shock into the island 
region, leading to higher particle and magnetic pressure and thus 
higher synchrotron emission which is more strongly boosted back 
into the minijet region co-moving frame C 2 . 

On the other hand, relativistic current sheets with no guide 
field have been found to develop a relativistic drift-kink instabil- 
ity (RDKI), which can disrupt the system before the part icles can 
be non- thermally accelerated (Zenitani & Hoshino 2008). To ex- 
plain TeV spectra in the flaring state of PKS 2155-304, a non- 
thermal power-law tail is needed in the electron distribution. A 
guide field has the effect of suppressing RDKI, allowing for ef- 
ficient particle acceleration. Models with a significant guide field 
(case II) can satisfy observational constraints, when radiation by an 
opposite minijet is taken into account. In fact, this effect is much 
more pronounced in case II, increasing the Compton dominance by 
2 orders of magnitude. Note, however, that these numerical stud- 
ies were done for pair plasma, while in our model electron-proton 
plasma is required. 

The amount of guide field in the minijets affects the spectrum 
in the soft X-ray band. This is independent of the slope of the non- 
thermal power-law tail (it is also true with no tail), but is related to 
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the ratio between synchrotron components produced in the minijet 
and island regions. In case I, the spectrum is hard, because emission 
from the island region is stronger due to stronger plasma compres- 
sion. In case II, the spectrum is soft, but still slightly harder than 
Chandra spectrum of PKS 2155-304. In the flaring state of this ob- 
ject, a harder-when-br ighter behaviour has be en observed in both 
X-ray and TeV bands JAharonian et alj|2009bh . This can be under- 
stood if the brighter flares are produced by the unguided minijets, 
while the fainter flares (and some part of the quiescent emission) 
come from the guided ones. 

An isolated event like a TeV flare in PKS 2155-304 should 
be associated with a significant, brief and temporary change in 
jet physical parameters. A single disturbance comoving with the 
bulk jet flow would cover a distance Aj- ~ r? ct c At ~ 
0.08(r jct /10) 2 (Ai/l d) pc. Thus, a ~ 4-day-long period of high 
activity would be related to a single global reversal of jet magnetic 
field travelling about 0.3 pc. If such a flare were triggered by a fac- 
tor external to the jet flow, it would produce a much longer activity 
period. 

An alternative scenario for the origin of minijets would in- 
volve kink instabilities developing in jets dominated by toroidal 
magnetic fields. During a flare, the jet would experience a brief 
global instability. A physical mechanism triggering it would be re- 
lated to some internal disturbance of the jet flow and not an external 
factor. Also, if the minijets arise from kink instabilities, the current 
sheet would have a more irregular structure and a case of two mini- 
jets aligned head-on is less likely. Thus, in this case it would be 
difficult to obtain a high Compton ratio, unless substantial external 
radiation fields are present in the reconnection region. 

Observations of PKS 2155-304 indicate that emission in the 
quiescent state may be of different origin than the bulk of emis- 
sion in the flaring state JAbramowski et alj2010h and that what we 
observe in the flaring state is a superposition of these two com- 
ponents JAharonian et alj|20 09b). If the quiescent state emission 
is of the same level during the flaring state, it might illuminate 
the minijets, providing some external radiation. Assume that the 
source of this emission is located downstream from the minijets 
zone and has the same bulk Lorentz factor Tjct = 10. The energy 
density of this external radiation of observed luminosity L b s in 02 
is u ext ~ £obsr2/(327rPj 1 ct r 2 c), where r is the distance between 
the stationary radiation source and the current sheet in the frame 
external to the jet. For an observed X -ray luminosity in the q uies- 
cent state L a \, s ~3x 10 45 erg s _1 JAharonian et al.ll2009ah and 
F2 = 10, we obtain u" xt ~ 10 rj~ 5 erg cm -3 . This is compara- 
ble with values shown in the upper left panel of Fig. [4] given that 
r ~ 10 15 cm, which is extremely close in terms of typical minijet 
sizeZ2 ~ 10 14 ~ 15 cm. It would be completely unimportant, if it ar- 
rived from a distance r ~ 10 18 cm, which is a typical location for 
the blazar zone in conventional models. Quiescent emission could 
be produced at very short distances from the cental black hole by 
minijets driven by kink instabilities. 

There are now indicatio ns that 7-ray emi s sion m ay be also 
anisotropic in bright blazars. ISavolainen et al.l ( 1201 Of) calculated 
comoving-frame viewing angles of a large sample of blazars 
(mostly Flat-Spectrum Radio Quasars), based on VLBA esti- 
mates of both Lorentz and Doppler factors. They have found 
that Fermi/LAT '-bright sources fell into a range of viewing angles 
40° < 9' < 110°, while those sources not detected in 3 months 
of Fermi/LAT monitoring had an almost uniform distribution in the 
full range of angles. This result cannot be easily explained without 
assuming internal anisotropy of the emitting region. The minijets 
model provides a physical structure that produces emission con- 



centrated in the direction perpendicular to the jet bulk flow, as mea- 
sured in the jet co-moving frame. On the other hand, variability 
constraints from Fermi/LAT are not strong enough to discriminate 
this from other possibilities. 



7 SUMMARY 

We have described a detailed model of minijets combining a dy- 
namical solution of relativistic magnetic reconnection with a weak 
guide field with calculations of non-thermal radiative processes in- 
cluding evolution of the electron distribution, pair-production opac- 
ity and Compton drag with an exact treatment of the Klein-Nishina 
cross section. Here are our main results: 

• TeV luminosities produced in model s using parame t er val - 
ues derived for a spherical blob model in iGiannios et alj (120090 
are much lower than those observed during fast TeV flares in PKS 
2155-304. This is because the minijet model presented here is lo- 
cated directly within the reconnection region of tiny opening angle. 

• Maximum 7-ray luminosity obtained for models with max- 
imum minijet region length allowed by observed variability 
timescale of 5 min is ~ 10 48 erg s _1 both for models with no or 
weak guide field and regardless of whether one takes into account 
radiation from an opposite minijet system. The luminosity can be 
limited either by opacity or strong Compton drag decelerating the 
minijet flow. External radiation from the opposite minijet allows 
one to obtain the maximum luminosity with higher Compton dom- 
inance. 

• SEDs matching the H.E.S.S. spectral fits for the flaring state 
of PKS 2155-304 differ widely with respect to the X-ray flux level. 
Simultaneous Chandra data are inconsistent with models including 
a significant guide field and excluding radiative interaction with 
an opposite minijet. In models with no guide field, radiation from 
misaligned minijets has to be taken into account in the X-ray band. 

• Our model can be fitted to PKS 2155-304 data only for jet 
magnetization ui > 25. Assuming bulk jet Loren tz factor Fj e t ~ 
10, th is is consistent with one-zone results of iBegelman et al.l 
l2008f) . The minijet concept allows one to reconcile modest bulk 
jet Lorentz factors with large local Lorentz factors in the TeV flar- 
ing region. 



APPENDIX A: SUMMARY OF RADIATIVE PROCESSES 

Synchrotron emissivity is calculated from the formula 
dCrusius & Schlic keiser 1986; Moderski et alj2003J) : 



3\/3 a T cu B I , 1 s-, I 
Jsyn(^) = -jz?^. / d-yn^jTZ 



(Al) 



47T 2 Vb J \3y 2 VB 

where Tl(x) = x 2 [K 1/3 (x)K 4/3 (x) - 0.6x(K 4/3 (x) 2 - 
K 1/3 (x) 2 )], u B = B 2 /(8tt) and v B = eB/(2-n-m e c). Since 
J 1Z(x)dx = 4v / 3~7r/81, one finds that the frequency-integrated 
formula is 



°TC , 2 \ 



(A2) 



The IC radiation emissivity is ca lculated from 
JAharonian & AtovarJll98ll ; lModerski et alj2005h : 



. , > 3 hcrr 

jic[y,n) = — e 

167T m e c 



1^4 



(7) 



f , u (e ) 
x / deo — /mO,6 m ) 



(A3) 
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where eo is the energy of the incident photon (in units of m e c 2 ), e 
is the energy of the scattered photon, «o(eo) is the energy density 
spectrum of incident radiation and /i is the cosine of the scattering 
angle, w — e/7, 6 M = 2eo7(l — A 4 )- F° r me / M function see Eq. 
JB2b . The frequency-integrated formula is: 



J'ic(m) = ^-(!-M) 2 / d7 7 2 "e(7) x 
x / deoWo(eo)/KN,^(b^) , 



(A4) 



where 



/kn, m (& m ) — 7J 



__5_ Jl (3 + 6,0(2 + 36,,) 
12 + 26„ 126 M (1 + 6 M ) 3 



+ (^-6)^ + 2) lri(l + ^ 



In the Thomson limit (6 M <C 1) /kn,/j(6 m ) ~ 1, thus 
Jic(m) . 3 2 uo_ 

JSYN 4 Mb 



(A5) 



(A6) 



When the incident radiation field is isotropic, we can average the 
above formulae over the scattering angle: 



3 h(T T 



jic(v) = -rr. e / dj 

16ty m e c J 7^ 



arc 



"e (7) 



deo 



Mo(eo) 



f(w, b) , (A7) 



j'ic = -5— / ^77 n e (7) / de Wo(eo)/KN(6). 



(A8) 



where 6 = 4eo7, 



/OM) = 1 + 



+ 



2tu 2 



2(1 -w) 6(1 -to) 6 2 (l-™) : 
6(1- to) ' 



+ 



2<c 



26(1 



/kn(6) 



6 + 



2 6(1-™; 
6 3 



In 



12(1 + 6) 



■ + 



6 3 



(^+6+^) 



(A9) 



(A10) 



In the Thomson limit, jic/jsyn ~ uo/ub- 

Radiative cool i ng o f electrons is calculated from 
JModerski et al.ll2003l.l2005l) : 



dx 



4<T T (7 2 - 1) 
3m e c 2 y/T 2 - 1 



UB 



deoMo(eo)/KN(6) 



(All) 



The pair-production absorption coefficient from a directed 
beam of ambi ent photons of energy ep = /i^o/(m e c 2 ) and energy 
density u$ is JGould & Schredenl 19670 



Kj-f{e, ji) — 



, Mo(eo) 



eo 



-ee 



(A12) 



where 



o 

v-,~,(x) = Jq°t (1-/3 2 ) x 



{Z-P A )^(j^)-W{2~?) 



(A13) 



and = sJT^TJx. 

© 0000 RAS, MNRAS 000, 000-000 



APPENDIX B: COMPTON DRAG IN KLEIN-NISHINA 
REGIME 

An integral formula for radiative force from the IC process tak- 
in g into account the Klein-Nishina cross-section was introduced 
in Blumenthal (1974) and used for studying Compton drag in 
iMadau & Th ompson (2000). We provide here an analytical for- 
mula, which is valid for an isotropic distribution of relativistic elec- 
trons of random Lorentz factor 7 2> 1, implying e^> eo, where eo 
is the energy of incident photons (in units of m e c 2 ) and e is the 
energy of scattered photons. I n this approximation, the di fferential 
IC cross-section is (E q. 20 in lAharonian & Atovanlll98ll : see also 
iModerski et~ai]|2005h : 



de dfl 



(eo,7) = 



3<tt 

167T7 2 eo 



/m( W A): 



(Bl) 



where e is the energy of scattered photons, /j, is the cosine of the 
scattering angle, w = e/7, 6 M = 2eo7(l — /tt), 



fn(w,bp) = 1 + 



■2w 



+ 



2w 2 



(B2) 



2(l-w) 6 M (1-™) 6 2 (l-w) 

The energy of scattered photons is limited to e < 76,,/ (1 + 6 M ). 

For a single scattering the electron loses energy of ~ em e c 2 
and gains momentum ~ —fiem e c in the direction of the incident 
photon. We calculate the average force per electron exerted by a 
directed photon beam of number density n p h(eo): 



= (7) 



de e 



deo n ph (e ) 
~ o°" T 7 2 / de u ph (e ) /kn,i(6) , 
where Uph(eo) = e m e c 2 n ph (e ), 



dfl /x ^ jn (eo,7) 



de dSl 



(B3) 



/kn,i(6) 



7 2 6 2 



de e 



22 + ^±^S-10 
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+2^— -Li 2 (-6) 



12(1 + 6) 2 



dw fn(w,b) = 
^Ml + 6)- 



(B4) 

= 76/(1 +b),b = 4e 7 and 



Mmax = 1 — iu/[2e 7(l - to)], e E 
Li 2 is the dilogarithm function. 

In the Thomson limit (6 <C 1), it can be shown that 
/kn.i(6) ~ 1 a nd so Eq. JB3b is consistent with Eq. (8) in 
ISikoraetal.l ( ll996h . 



APPENDIX C: RADIATION ENERGY DENSITY FROM A 
STATIONARY PATTERN OF RELATIVISTIC FLOW 

A recurring problem in this work is to calculate the energy den- 
sity of synchrotron and external radiation in a frame co-moving 
with a relativistic flow. The emitting region is stationary in frame 
O and the emitting fluid is characterized by bulk Lorentz factor 
r em and co-moving emissivity j'{v' em ). Consider a small element 
of the emitting region of volume AV, observed from distance r, 
with fluid velocity making angle 8 with the direction away from the 
observer. At the observer's position, there is a flow of Lorentz fac- 
tor r o bs and velocity parallel to the velocity of the emitting fluid. 
In O, the radiation energy density is Au(v) — j(v)AV/(r 2 c). 
Emissivity in O is given by j(u) = / D 2 m j'(V~ I lv), where 
2?om = [r cm (l + /3 om cos^)]^ 1 . A transformation of energy 
density to the fluid frame at observer's position is Au'(v') — 
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T>' ohs Au(D' obs V), where V' ohs = T obs (l + /3 obs cosfl). In ef- 
fect, we have obtained: 



AuV) = v oha n 



f (k 



^cm^' J /\y 



(CI) 



c r^ 

In application to radiation energy density in the minijet co- 
moving frame, we would take r o b s = T2 and r em = Tg for ra- 
diation emitted in the island region or r cm = T2 for synchrotron 
radiation emitted in the minijet region. In the case of r cm = r o b s , 
formula dClb reduces to 



Au'(i/) = V c 



f («/') AV 



(C2) 
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